Critical scaling of thermal transport in model A dynamics for a superconductor 
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We consider the scaling of thermal transport in the presence of electric and magnetic fields near 
the finite temperature transition from the metallic normal state to the superconductor. We do so 
with fully relaxational, model A, dynamics for the order parameter and particle-hole symmetry. 
This enables us to determine the exact scaling dimension of the heat current operator and hence 
critical exponents for all transport coefficients in terms of the (approximately) known values of the 
correlation length exponent and the dynamic scaling exponent. In particular, we determine the 
critical behavior for the Nernst coefficient. 



I. INTRODUCTION 

The discovery of high-temperature superconductors 
has led to a renewed interest in superconducting fluctu- 
ations. In these materials, the effect of fluctuations may 
be strong, with a relatively large critical regime, due to 
their short coherence lengths, quasi- two-dimensionality 
and high transition temperatures Q. Neglecting fluc- 
tuations of the gauge potential, the transition to the 
superconducting state is described by a complex order 
parameter and is of the XY universality class (in three 
dimensions, due to interlayer coupling). This description 
holds well in strongly type-II superconductors except ex- 
tremely close to the critical temperature T c where one 
crosses over to the critical behavior of a charged super- 
conductor 0| . The description of dynamic properties, 
such as transport, requires coupling this thermodynamic 
description to a dynamic equation, which for supercon- 
ductors is usually assumed to be model A in the classifi- 
cation of Hohenberg and Halperin Q . 

Of particular interest in this context are the re- 
cent measurements of the Nernst coefficient in the 
cuprates Q|. The large Nernst signal observed in these 
experiments above T c is likely the contribution of su- 
perconducting fluctuations. This interpretation provides 
a quantitative agreement with experiment in overdoped 
samples in the regime where a Gaussian approximation is 
applicable [|| , as well as in the vortex liquid regime |(| . In 
underdoped samples it requires a stronger effect of fluc- 
tuations, which would have important implications for 
the physics of the pseudogap regime pfl . 

In this paper, we focus our attention on the scaling 
behavior of thermal transport in the XY critical regime. 
The transport coefficient of interest for the Nernst effect 
is the transverse thermoelectric response a xy p|. For 
this, we consider the scaling of the heat current operator, 
which also allows us to treat the thermal conductivity on 
the same footing. As with the conductivity 0, @ j it is 
necessary to specify the dynamics associated with the or- 
der parameter for analyzing thermal transport. We then 
find the scaling dimension of the heat current operator, 
and obtain general scaling forms for thermal transport. 
The critical behavior of both the thermal and thermo- 
electric transport coefficients is then deduced. 



We follow the usual assumption for a superconductor, 
and consider relaxational dynamics for a non-conserved 
order parameter, or model A dynamics 0,0. We fur- 
ther assume that the model is particle-hole symmetric 
which allows much progress to be made. When needed, 
the order parameter may be coupled to an electromag- 
netic field, whose dynamics we ignore, keeping our calcu- 
lation in the XY regime [Til ]. 

The main result of this paper is simply stated: the heat 
current operator in model A, 



dvb* 



(i) 



does not acquire an anomalous dimension. This result is 
surprising on two counts: First, the heat current is not 
a conserved quantity in model A and second, the scal- 
ing dimension of the heat current does not involve the 
dynamic critical exponent z, despite its dynamic nature. 
Technically we will find that there is a "hidden" conser- 
vation law at play in the model. We note that a similar 
absence of an anomalous dimension is well known for the 
electric current Q, 



j cx —iijj*\/tjj + c.c 



(2) 



In this case, however, there are conserved (solenoidal) 
supercurrents in the purely static theory and thus the 
result is less mysterious. 

With this result, critical exponents of various thermal 
transport coefficients may now be deduced. In partic- 
ular, we find that to linear order in magnetic field, the 
transverse thermoelectric response a xy scales in three di- 
mensions as a xy /B cx , where £ is the (temperature- 
dependent) correlation length of the superconductor. 
Whether this prediction may be observable in the critical 
3D-XY regime of cuprates remains to be seen. 

Below we present our calculations in detail: After 
defining the model in Sec. [HI we obtain the anomalous 
dimension of the heat current operator (or rather, its ab- 
sence thereof) in Sec. lIIII We briefly address the problem 
of introducing a temperature gradient in Sec. IIVI The 
scaling relations for thermal transport and their conse- 
quences are discussed in Sec. E] which is followed by a 
brief summary. 
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II. MODEL A 

In this paper, we consider model A dynamics for a 
complex superconducting order parameter. The order 
parameter dynamics in this model is given by the stochas- 
tic equation 



84, _ SF 
with the Landau-Ginzburg- Wilson free energy 



F 



dx(r |Vf +u |V| 4 + |VV| 5 



(3) 



(4) 



In Eq. Pjl. r is the relaxation rate for the order param- 
eter, which we assume in this paper to be real (making 
the model particle-hole symmetric). Thermal fluctua- 
tions are introduced with £, which is a Gaussian white 
noise with correlator 



<C* (x, t)C(x', t')) = 2r T<5(x - x')5(t - t') 



(5) 



Unlike the usual convention in critical dynamics, we have 
explicitly expressed the temperature T, which will be use- 
ful below. The correlation function of the noise ensures 
that the probability that an order parameter configura- 
tion ip will occur is P{ip} = Z^ 1 cxp[— F{-ip}/T] (where 
Z = fD{i,}exp[-F{^}/T}). 

The main object we consider within this model is the 
operator of heat current (Q. This form may be under- 
stood in terms of the transport of free energy. In Sec. II VI 
we note how it arises from considerations involving the 
application of a temperature gradient to the system. It 
also arises in derivations of time-dependent Ginzburg- 
Landau equations from microscopies. Note that the heat 
current operator is a dynamic object involving a time 
derivative (unlike the electric current); thus its proper- 
ties are dependent on the choice of dynamics. 

When particle-hole symmetry is broken, the relaxation 
rate Tq becomes complex, and model A admits non- 
dissipative, or reactive, terms. While we expect such 
terms not to contribute to heat transport, the question 
of incorporating this expectation into the critical dynam- 
ics treatment of the heat current is not fully resolved. 
In this paper, we therefore proceed with the assumption 
particle-hole symmetry. 



III. SCALING DIMENSION OF j Q 

In this section, we consider the scaling dimension <1q of 
the heat current operator. The dimension cLq will then 
allow us to write scaling forms for properties involving 
the heat current, as we do in Sec. for thermal and 
thermoelectric transport. Naively, one may expect the 
emergence of an anomalous dimension for this operator 
(i.e., a deviation from its Gaussian value). In particular, 
because of the explicit time derivative in [Eq. 



one may expect the anomalous dimension to involve the 
dynamical critical exponent z. In this section we show 
that these expectations are false; we find dg = d + 1 in d 
dimensions, with no anomalous dimension. We do so first 
by an exact argument, then by sketching the calculation 
of cLq to second order in an e expansion. 

We begin by replacing the dynamic stochastic equa- 
tion with an effective action for calculating dynamic cor- 
relators |T^ . Consider the Langevin equation given by 
Eq. © (in this section we set T = 1). The noise term 
has a Gaussian distribution, 



P{((x, t)} cx exp 



2r 



(6) 



Enforcing Eq. @ by a S function, expectation values, 
when averaged over the noise distribution, take the form 



•) c cx J DijjDC 



00 

The integral over the noise may now be performed, yield- 
ing an effective action for ip, 



where 



D4> 



/ d d xdt 

2r 



d± 5F_ 

dt °5ip* 



(8) 



(9) 
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The last step uses the fact that the cross term in the ac- 
tion is a full derivative and hence vanishes upon integra- 
tion. In the above, we have ignored the Jacobian term, 
which is irrelevant for the renormalization group analy- 
sis |l2j. The problem of calculating dynamical correla- 
tion function is thus expressed in terms of a new effective 
action S e g, which may be viewed either as a quantum me- 
chanical action (the view adopted here), or equivalently 
as a classical statistical mechanical action in d+ 1 dimen- 
sions. 

The exact argument regarding the dimension of the 
heat current operator is based on the following obser- 
vation: The heat current, Eq. is proportional to the 
conserved momentum density of the effective Lagrangian, 



d(d ip) T \ dt 



(10) 



This observation essentially sets the dimension of the 
heat current operator; it does not have an anomalous 
dimension because it is conserved. A similar argument, 
but in the purely static theory, fixes the dimension of the 
electric current operator as d e = d — 1 . 

A standard argument turns the conservation law into a 
Ward identity. Consider a time-ordered correlation func- 
tion of Toi with a set of operators ip (whose dimension is 
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known) , 



(11) 



(Here, we use x to denote both x and t in the arguments 
of the fields.) Using d^T^ = (as T Mi is a conserved 
current for each of its components i), we only need to 
account for the time-derivative on the step functions due 
to time-ordering, 

^<TT^(x)V>(zi)...#r„)) (12) 

n 

= J2 5 (t- **) Wsi) ■ ■ ■ PM, ^)] . . . V(*„)> 

3=1 

n 

= -*£ ^ (d+1) (^ - ax)M*i) • ■ ■ dMxj) . . . if,(x n )). 

3=1 

The dimension of the heat current operator, c?q = d+ 1, 
may be deduced from Eq. I|12[l by comparing the term 
with (i = in the first line (involving the required heat 
current operator) with the last line of the equation (where 
all the dimensionalities are known). The dimensionality 
of the fields ip(xi) cancel between the two sides, as does 
the factor z between the time derivative in the first line 
and the time component of the S function in the last line. 
The spatial part of the S function and the additional 
spatial derivative in the last line give the dimension of 
the heat current operator. 

As the result above is not very transparent in terms of 
its physical content it is instructive to see how this works 
in an expansion about the Gaussian fixed point; indeed, 
this is how we came across this result in the first place. 
Specifically we consider the expansion in e = 4 — d to sec- 
ond order (the first non-trivial order) for the dimension 
of the heat current. We find that the coefficients of this 
expansion vanish, thus corroborating the exact argument 
above. 

The calculation proceeds by adapting the method for 
computing dimensions of composite operators outlined 
by Wilson and Kogut 01 to the dynamic theory. For 
convenience, the complex order parameter ip is repre- 
sented in this calculation in terms of its real and imag- 
inary components, ip — <j)\ + ifo, and the heat current 
operator is then 



(13) 



We consider the following correlation function involving 
the heat current operator 

Oh.)./^**-*^ (14 , 

x (0 2 (x,t)^(x',t')j Q (O,O)). 

Following Wilson and Kogut, the correlation function has 
been chosen both for its simplicity and because it avoids 
the mixing in of operators involving total derivatives. 



By simple counting of dimensions, the correlation func- 
tion may be expressed (in terms of an unknown scaling 
function F) as 



2d+2z~2d s 



(15) 



Here, £ is the correlation length, cZq is the unknown di- 
mension of the heat current operator, and d s = (d — 2 + 
r])/2 is the dimension of the operator This result 
may be expressed in terms of the {<j>i<j>i) static correlator 
at q = 0, 7 ,_1 cx £ 2 ~''. Expanding to linear order in q 
and u), we thus have 



Uocqwr( d Q+"- d - 3 - 3z )/( 2 - 1 '). 



(16) 



This result will be compared with a perturbative calcu- 
lation of U to obtain the result for dQ. 

We briefly outline the calculation of U. The pertur- 
bative method for the dynamic case (see, e.g., Ref. prf) 
begins by presenting the fields as 



&(k,w) =^ (0) (k,c)-uGo(k,c)]r J 



dk' dk" du' duj" 



(27r) 2 ( d 



+i) 



x 0j(k', u')<t>j(k", w")<Mk - k' - k", u - <J - u"). (17) 



Here, 



Go(k, w)£j(k, u)/r is the free field, and 
G (k,w) 



r + k 2 — iuj/T 



(18) 



is the free propagator. Equation (|17fl is then used to 
generate the different diagrammatic contributions to the 
desired order in u. In contrast with the static perturba- 
tion theory, such diagrams contain lines corresponding to 
both the propagator Go and the correlator 

(#(k, W )0f (k>'))= (19) 
(2n) d+1 S{uj + uj')S(k + k')<%G (k, u), 

where Go(k, ui) = (2/w)ImGo(k, ui). The diagrams which 
arise to second order in u are presented in Fig. ^ 

At Gaussian order [see Fig. [Ja)] and to linear order in 
q and lo, we have 



U 



4qw 



(20) 



To first order in u we have the diagram in Fig. ^b), 
but its contribution vanishes due to symmetry consider- 
ations [as does the second order diagram in Fig. ^c)]. 
(Diagrams in which a closed loop appears on one of the 
propagators are included by the renormalization of r.) 
To second order in u, we have the diagrams in Fig. ^c)- 
(c). Of these, we have illustrated in Fig. ^d) the fact 
that each of these diagrams is in fact a set of diagrams 
with different placements of propagator and correlator 
lines. Note that the external legs in the different dia- 
grams may be either propagator or correlator lines. For 
this reason, in contrast with the static case, we calculate 
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(a) , (b) , (c) 




FIG. 1: Diagrams appearing in the correlator U to second 
order in u. The heat current vertex j*^ is denoted by a black 
circle. A line with an arrow denotes the bare propagator Go, 
and a line with an empty circle denotes the bare correlator 
Co [in (b), (c), and (e) the diagrams are shown without pre- 
senting all possible arrangements of propagator and correlator 
lines]. Diagram (a) is the bare (or zeroth order) contribution. 
Diagrams (b) and (c) are first and second order diagrams, 
respectively, which have a vanishing contribution. Diagrams 
which contribute to second order are given in (d) and (e). 



all diagrams appearing to a given order, not only those 
which are irreducible with respect to the external legs. 

In the calculation of the diagrams, we are interested in 
the infrared divergence when r — > 0. The diagrams are 
also ultraviolet divergent, which we re gula rize by dimen- 
sional regularization (see., e.g., Ref. [l5j|l. The actual 
calculation of these diagrams is rather lengthy, and we 
will spare the reader the details. Very briefly, the various 
diagrams are evaluated using the value of u at the fixed 
point, which to order e is u* — tt 2 e/5. It is convenient to 
express the result in terms of the exponent r\ (to second 
order in e, r\ = e 2 /50) and the coefficient c = 61n| — 1, 
which appears in the result for the dynamical critical ex- 
ponent, z — 2 + cq jlfij |. The result is 



U = U (l--(c+l)rilnA oc r- 4 - 3 ( c+1 )"/ 2 . (21) 



A comparison of Eqs. (|16|l and l|21|l gives (1q = d + 1 to 
second order in e. We thus do not find an anomalous di- 
mension for the heat current operator in the e expansion, 
in agreement with the general argument given above. 



IV. INTRODUCING A TEMPERATURE 
GRADIENT 

In previous sections, the temperature T was constant 
throughout the sample. A constant temperature may 
be simply absorbed by an appropriate rescaling of vari- 
ables, so that it does not appear explicitly in the problem. 
Indeed, by rescaling the order parameter, free energy, 
and noise using tp — ► yTijj, F — ► TF, and C - * C/VT, 
Eqs. m a<y be rewritten without an explicit tem- 

perature (this also requires rescaling the coefficient of the 
quartic term in the free energy). This is the form tradi- 
tionally used in the literature on critical dynamics 

For thermal and thermoelectric transport, however, we 
should also consider the application of temperature gradi- 
ents in the sample. For this purpose, we have retained the 
explicit temperature dependence of the noise correlator, 
Eq. JSJ). Other parameters of the problem may also have 
an implicit temperature dependence; however, at least 
for linear response to a temperature gradient, the de- 
pendence of parameters on temperature is not important 
for transport and is subsequently ignored. The reason 
for this is that if we ignore the temperature dependence 
in the noise, the system then has spatially-dependent 
local couplings (through their temperature dependence) 
but remains in isothermal equilibrium; thus it carries no 
transport currents. 

Rescaling Eq. ©, as discussed above, leads to the fol- 
lowing equation of motion, to linear order in the temper- 
ature gradient, 



- r [r ^ + 2u |</>|V - V 2 V - (WO • (VT)/T] + c 

The dependence of the correlator (J5J on temperature has 
been exchanged for a term proportional to VT in the 
stochastic equation. The equation may still be recast in 
the form Q, provided an additional term, of the form 
— j*2 • A, is added to (the integrand of) the free energy, 
Eq. 0. Here, is the heat current, Eq. (Q. It is cou- 
pled to a "gravitational" vector potential A. The notion 
of a gravitational field was introduced by Luttinger \vfy 
as a fictitious device for studying temperature gradients 
(see also Ref. Here, we briefly discuss this point, 

noting the analogy between the treatment of electric and 
thermal currents. In particular, Einstein relations for 
thermal currents relate the response to (VT)/T with the 
response to a gravitational force E — — V</>, where <f> is 
Luttinger's gravitational field. By performing a gravi- 
tational gauge transformation, the gravitational force is 
expressed in terms of a gravitational vector potential, 
E = —dA/dt. To linear order in the gravitational force, 
the gravitational vector potential A is then coupled to 
y« in the free energy of model A. The equation of mo- 
tion 0, coupled with the Einstein relations, reproduces 
Eq. (22). 

This understanding of how a thermal gradient enters 
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the free energy is useful on several fronts. First, it verifies 
the form of the heat current operator, Eq. (JIJ . Second, it 
replaces temperature gradients with gravitational fields. 
Luttinger's original motivation was to express thermal 
transport coefficients in terms of the corresponding cor- 
relation functions. In the next section, we follow a similar 
strategy for writing scaling forms for the thermal current. 

V. THERMAL TRANSPORT 

In this section, we consider the critical behavior of 
thermal transport coefficients in model A dynamics. Us- 
ing the result for the scaling dimension of the heat current 
from Sec. IIIII dq = d + 1 , we write the general scaling 
relations for the electric and heat currents, in presence of 
an electric, magnetic, and gravitational fields, 



-d+1. 



-1+z 



+z 



(23) 
(24) 



These equations are extensions of known results for elec- 
tric transport (see, e.g., Ref. Q). 

We will use these scaling relations here to deduce the 
critical behavior of various linear response coefficients. 
We note that the assumption of particle-hole symme- 
try implies that a xy = a xx — n xy — (it is necessary 
to break particle-hole symmetry to discuss the critical 
properties of these transport coefficients). For the con- 
ductivity, as well as for the magnetic susceptibility, the 
known results jlj may be confirmed, 



z — d 



x « 6 



4-d 



(25) 



Our result for the heat current operator now enables us 
to make analogous statements for transport coefficients 
involving the heat current operator, a xy and n xx . 

We begin with the transverse thermoelectric response 
a xy , which is the coefficient of interest in studying the 
Nernst effect. More precisely, we will consider the re- 
sponse to linear order in the magnetic field. This may be 
calculated either from the heat current response to an ap- 
plied electric and magnetic field, or as the electric current 
response to the applied gravitational and magnetic fields. 
In both cases, the current which is obtained is the total 
current, from which magnetization contributions, which 
do not contribute to transport, must be subtracted |18| . 
For the total currents, we find 



Jtot Jtot s2+z-d 

EB (V0)5 



(26) 



The magnetization component of these total currents is 
given by 0] 

i Q = M x E 



Jn 



= M X V(j>. 



(27) 



Using the result for the magnetic susceptibility above, 
we see that the situation is different dependent upon the 
value of the critical dynamical exponent. 



For d > 4 mean-field exponents hold, and the dynami- 
cal critical exponent is z = 2. In this case, the total and 
magnetization pieces of the current diverge with identi- 
cal power laws (cf. Ref. p|). The difference of two terms 
with identical singularities may be of the same singular 
behavior or a weaker one. Thus, the critical dynamics 
analysis presented here does not determine conclusively 
the power of a xy when z = 2. For Gaussian fluctuations, 
a xy does have the same singularity as the total currents 
and magnetization contributions, a xy /B oc £ 4_d jjj. It is 
quite possible that the same occurs also in two dimension 
(where z = 2 as well), suggesting that a xy /B oc £ 2 above 
the Kosterlitz-Thouless transition temperature. How- 
ever, a separate analysis is required for this case which 
we defer to a future work. 

For three dimensions, which corresponds to the actual 
transition in the cuprates due to interlayer coupling, the 
analysis of Ref. 16J gives z = 2 + crj > 2. In this case, the 
magnetization currents are less singular than the total 
currents from which they are subtracted, and hence the 
result for a xy to linear order in B, 



a xy /B oc £ 



2+z- 



(28) 



becomes unproblcmatic .19]. We note that the dimcn- 
sionless ratio Ta xy /M (where M = \B is the magneti- 
zation) diverges at the critical point as £ z ~ 2 . 

Finally, we consider the thermal conductivity, for 
which we find 



-2+z-d 



(29) 



In particular, the thermal conductivity will have a singu- 
lar, but non-diverging contribution in model A in three 
dimensions. In contrast, Vishveshwara and Fisher argued 
in a recent paper M that the thermal conductivity is an- 
alytic at the critical point using a model C formulation 
(which involves a conserved energy density in addition 
to a non-conserved order parameter). There are several 
things to note about this difference in the critical behav- 
ior obtained in the two works. First, the coupling of the 
energy variable to the order parameter is known to be 
irrelevant (as a < with n = 2 components of the order 
parameter), and we expect the critical behavior of model 
C to be equivalent to that of model A |2l|. Second, we 
note that singular behavior is already obtained for the 
thermal conductivity for Gaussian fluctuations in model 
A p|, a result which may also be obtained from micro- 
scopies (which are energy conserving). Finally, it seems 
to us that in Ref. [2(j , a treatment of the order parameter 
contribution to the heat current is absent. 



VI. FINAL COMMENTS 

Our results for the critical behavior of the transverse 
thermoelectric response a xy and the thermal conductiv- 
ity k X x are based on the observation that the heat current 
operator has no anomalous dimension. It is worth noting 
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that both electric and heat currents share this property; 
in particular they are independent of the dynamical crit- 
ical exponent z. Readers familiar with the lore of quan- 
tum critical phenomena should note that conserved cur- 
rents do exhibit z in their scaling dimensions at critical 
points; it is their corresponding densities which do not. 
Thus the meaning of the phrase "absence of anomalous 
dimension" is different in these two cases. 

We briefly comment on the possibility of observing our 
results in experiment. In the moderately two dimensional 
cuprates, such as YBCO, there is evidence for a rela- 
tively large 3D-XY critical regime. In this regime our 
results for a xy , to linear order in the field, at tempera- 
tures approaching T c can be tested via measurements of 
the Nernst effect and the conductivity. Another possibil- 
ity, instead of attempting to extract critical exponents, 
is to consider the dimensionless ratio Ta xy /M. This ra- 
tio equals 1/2 in the region of Gaussian fluctuations [jj. 
As the temperature is decreased towards T c , the value 
of this ratio increases, eventually diverging as £ z ~ 2 . In 
the highly two dimensional cuprates such as BSSCO the 
scaling theory is on somewhat weaker grounds due to the 
degeneracy between the dimensions of the total current 



and the magnetization. However, absent an exact cancel- 
lation, it should still correctly predict that a xy to linear 
order in B diverges with the same critical exponent (i.e. 
~ £ 2 ) as the magnetization as the Kosterlitz-Thouless 
transition is approached (except perhaps for a confluent 
logarithm at the lower critical dimension.) 

In summary, we considered the critical scaling of ther- 
mal transport near the finite temperature transition to 
the superconducting state, within the framework of crit- 
ical model A dynamics. This leads to specific predic- 
tions for the critical exponents of thermal transport co- 
efficients, and in particular for the Nernst effect. 
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